Fluctuating spin g-tensor in small metal grains 
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In the presence of spin-orbit scattering, the splitting of an energy level e M in a generic small metal 
grain due to the Zeeman coupling to a magnetic field B depends on the direction of B, as a result of 
mesoscopic fluctuations. The anisotropy is described by the eigenvalues gr| (j = 1, 2, 3) of a tensor 
Q, corresponding to the (squares of) (/-factors along three principal axes. We consider the statistical 
distribution of Q and find that the anisotropy is enhanced by eigenvalue repulsion between the gj . 

PACS numbers: 71.24.+q, 71.70.Ej 



With the advance of nanoparticle technology, it has 
become possible to resolve individual energy levels for 
electrons in ultrasmall metal grains. Recent experiments 
addressed their Zeeman splitting under the application 
of a magnetic field B [jf] || . The splitting of a level e M 
is described by a (/-factor, fc M = ±i/isgi? z , where [is is 
the Bohr magneton. A free electron has g — 2, but in 
small metal grains the effective (/-factor may be reduced 
as a result of spin-orbit scattering [Q . In order to study 
this reduction, Salinas et al. || have doped Al grains 
(which do not have significant spin-orbit scattering) with 
Au (which has). For small concentrations of Au, the 
effective (/-factor was seen to drop from 2 to around 0.7. 
Even lower values g ~ 0.3 were reported in experiments 
on Au grains (|]. 

For disordered systems with spin-orbit scattering, the 
splitting of a level does not only depend on the mag- 
nitude of the magnetic field B, but also on its direction. 
Hence, an analysis in terms of a "(/-tensor" is more ap- 
propriate ||. To be precise, the Zeeman field splits the 
Kramers' doublet — > ± fe M with 

Sel = (ii B /2) 2 B ■g iM -B, (1) 

where £/ M is a 3 x 3 tensor. In the absence of spin-orbit 
scattering, the tensor is isotropic, {Q^)ij = 4(5y . The 
effect of spin-orbit scattering on is threefold: It leads 
to a decrease of the typical magnitude of Q^, it makes 
the tensor structure of important (i.e., it introduces 
an anisotropic response to the magnetic field B), and it 
causes to be different for each level e M . Hence <? M 
becomes a fluctuating quantity, and it is important to 
know its statistical distribution. The latter problem was 
addressed in a recent paper by Matveev et al. ||, how- 
ever without considering the tensor structure of . The 
anisotropy of the g-tensor is a measurable quantity and 
we here consider the distribution of the entire tensor Q^. 
The distribution P(G^) is defined with respect to an en- 
semble of small metal grains of roughly equal size. The 
same distribution applies to the fluctuations of as a 
function of the level e M in the same grain. 

In general, has a contribution G spin from the mag- 



netic moment of electron spins, and a contribution Q° 
for the orbital angular moment of the state |Vv)- In Ref. 
H , the typical sizes of both contributions were estimated 
as Q spin ~ t so A and G orh ~ i/L, where r so is the mean 
spin-orbit scattering time, L is the grain size, A cx L~ 3 
is the mean level spacing, and £ -C L is the elastic mean 
free path. We restrict ourselves to the spin contribution 
(z spm , which should be dominant for small grain sizes ||, 
provided t so does not depend on system size, as should 
be the case for the experiments of Ref. 0. When or- 
bital contributions are important, the anisotropy of Q 
will be affected by the shape of the grain. In that case, 
our main conclusions apply only to a roughly spherical 
grain. As the typical magnitude of Q (we drop the super- 
script "spin" and the subscript [i if there is no ambigu- 
ity) depends on the microscopic parameters t so and A, 
which are in most cases not known accurately, we choose 
to have the typical magnitude of Q serve as an external 
parameter in our theory 

We first present our main results. With a suitable 
choice of the coordinate axes ("principal axes"), the ten- 
sor Q can be diagonalized. Writing its eigenvalues as g 2 
and denoting the components of the magnetic field along 
the principal axes by Bj, j = 1, 2, 3, Eq. ([!]) takes a par- 
ticularly simple form, 

K = idlBl + g\B\ + g\Bl). (2) 

We refer to the numbers g\, g2, and g^ as principal g- 
factors. For a generic metal grain of a cubic material, 
rotational symmetry implies that, for a given level e^, 
the positioning of the principal axes is entirely random 
in space, as long as they are mutually orthogonal. Hence, 
it remains to study the distribution ^(31,32,33) of the 
principal g-factors g\, 32, and 173 for the level e^. Our 
main result is, that for sufficiently strong spin-orbit scat- 
tering, P{gi, <?2, 33) is given by the distribution 

p(gi, 92, 33) « n \ g f -sun e_39 ' /2</> ' ( 3 ) 

i<j i 

where g 2 = \{g\ +g% +93) is the average of (2fe M / ' ^jlbB) 2 
over all directions of B and (g 2 ) is its average over the 
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FIG. 1. Average of the squares of principal g- factors ver- 
sus spin-orbit scattering strength A, obtained from numerical 
simulation of the random matrix model (^) with N = 100. In- 
set: g\, Q2, and gs for a specific realization. We have included 
the sign of gi; see the discussion below Eq. (Js|). 

ensemble of grains. In random matrix theory [Q, this 
distribution is known as the Laguerre ensemble. With- 
out loss of generality we may assume that g\ < g\ < g\. 
Figure [l] shows the averages and a realization of the 
principal (/-factors g±, g%, and g$ for a specific sample, as 
a function of a parameter A ~ (tsqA)" 1 / 2 measuring the 
strength of the spin-orbit scattering. (A formal defini- 
tion of A in a random-matrix model will be given below.) 
From the figure, one readily observes that, typically, the 
three principal (/-factors differ by a factor 2-3. This im- 
plies that, in spite of the average rotational symmetry of 
the grains, the response of a given level to an applied 
magnetic field is highly anisotropic because of mesoscopic 
fluctuations. The mathematical origin of this effect is the 
"level repulsion" factor \gf — g?\ in the probability distri- 
bution (^), which signifies that, to a certain extent, 
can be viewed "random matrix" . 

Let us now turn to a more detailed discussion of our 
results. Without magnetic field, the Hamiltonian 7i 
of the grain is invariant under time-reversal, so that 
all eigenstates come in doublets IVv) and \Tip^}, where 
Tip = i<j2'<p* is the time-reversal operator. To study the 
splitting of the doublets by a magnetic field, we add a 
term [IbB ■ a to Ji, a — (o"i, 02, 0-3) being the vector of 
Pauli matrices. From degenerate perturbation theory we 
find that a level e M is split into £ M ± <5e M , with fe M of the 
form (|l|) . For the real symmetric 3x3 matrix one has 

Q„ = G£G M , (4) 
where G M is a real 3x3 matrix with elements 

+i(G^) 2i = -2(T^h-|Vv) 

(G„) 3J = 2<</vh# M >, (5) 



We use random-matrix theory (RMT) to compute the 
distribution of Q^. In RMT, the microscopic Hamiltonian 
Ti is replaced by a 2N x 2A~ random hermitian matrix H, 
where at the end of the calculation the limit N — > 00 is 
taken. (The factor 2 accounts for spin.) The wavefunc- 
tion ipf_i(r) is replaced by an A^-component spinor eigen- 
vector ippn of H, where n is a vector index. To study the 
effect of spin-orbit scattering, we take H of the form 



if (A) = s® a 2 +i- 
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(6a) 



where S (Aj ) is a real symmetric (antisymmetric) N x N 
matrix with the Gaussian distribution 

P(S)oce^ 2 / 4ArA2 > trSTs , (6b) 

P(A ) OC e -(- 2 /4A^)trAjA 3; j = 1)2 3 

The Hamiltonian H(X) is similar to the Pandey-Mehta 
Hamiltonian used to describe the effect of time-reversal 
symmetry breaking in a system of spinless particles 
||. In Eq. (6b), A is the average spacing between the 
Kramers doublets near e = 0. The amount of spin-orbit 
scattering is measured by the parameter A ~ (r so A) -1 / 2 
. The case A = corresponds to the absence of spin- 
orbit scattering, when H = S is a member of the Gaus- 
sian Orthogonal Ensemble (GOE) of random matrix the- 
ory. The case A = (AN) 1 / 2 corresponds to the case of 
strong spin-orbit scattering, when H is a member of the 
Gaussian Symplectic Ensemble (GSE). The ensemble of 
Hamiltonians H(X) corresponds to a crossover from the 
GOE to the GSE. Similar crossovers were studied previ- 
ously in the literature, in particular for the cases GOE- 
GUE and GSE-GUE (GUE is Gaussian Unitary Ensem- 
ble) 

The distribution of the tensor for an eigenvalue e M of 
the matrix H(X) is related to the statistics of eigenvectors 
of H(X) in this crossover ensemble. To deal with the 
twofold degeneracy of the eigenvalue e M , we combine the 
two A^-component spinor eigenvectors ip^ and Tip^ into 
a single A-component vector of quaternions ip = (ip, Tip) 
PlvM- The quaternion vector ip can be parameterized as, 
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(7) 



where the Uk are quaternion numbers with tvu\ui = 2Ski 
("quaternion phase factors"), the (pk are A^-component 
real orthonormal vectors, and the are positive num- 
bers such that J2k a 1 = 1 (fe, Z = 0,1,2,3). A eigen- 
vector in the GOE corresponds to ctQ = 1, a.\ = a.2 = 
as = 0, while an eigenvector in the GSE has typically 
ccq w a.1 ~ ^2 ~ a 3 Ri -i. A similar parameterization has 
been applied to the GOE-GUE crossover || . Orthogonal 
invariance of the distributions of S and Aj , together with 
the freedom to choose the overall quaternion phase of ip, 
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give a distribution of the Uk and 4>k that is as random 
as possible, provided the above mentioned orthogonality 
constraints are obeyed. Hence, all nontrivial information 
about the eigenvector statistics is encoded in the num- 
bers ctk- Substitution of the parameterization (ffl) into 
Eq. (§) yields 



0/2,2 2 2\ 

r,/ 2 2,2 2\ 

92 = 2(a -a 1 +a 2 - a 3 ), 

r,/ 2 2 2 i 2 \ 

33 = 2(a -a 1 -a 2 + a 3 ). 



(8) 



While the squares a 2 (k = 0, 1, 2, 3) are all positive, the 
principal 3-factors as given by Eq. (||) can also be nega- 
tive. Permutations of the at alter the signs of the indi- 
vidual gj, but not of their product 310233- [The product 
313233 = det G also follows from Eq. (g); one verifies that 
it does not change when \ip) is replaced by a linear com- 
bination of \tp) and Without loss of generality, 



3§, and that 32 and 33 
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we may assume that gf < g\ < 

are positive. Then equation (||) provides the constraint 
32+33 < 2+31, which poses a bound on the occurrence of 
negative values for the product 313233- We conclude that 
all information on the eigenvector statistics in the GOE 
GSE crossover is encoded in the magnitudes of 31, 32, 
and 33 and the sign of their product. Since for the level 
splitting 5e^{B) only the squares g 2 are of relevance, we 
disregard the sign of 313233 in the remainder of the pa- 
per. The sign of 313233 may be determined in principle, 
however, by a spin-resonance experiment ]Ti| ]. 

In order to calculate the distribution P (31,32,33) one 
has, in principle, to carry out the same program as was 
done in Refs. (jlO 11 for the GOE-GUE crossover. How- 
ever, it turns out that in the present case the calcula- 
tion is considerably more complicated. This can already 
be seen from the mere observation that the wavefunction 
statistics in the GOE-GSE crossover is governed by three 
variables 31, 32, and 33, whereas in the case of the GOE- 
GUE crossover only one variable was needed |jl^-|l2| . In 
the field-theoretic language of Ref. |ll[] , one has to use a 
nonlinear sigma model of 16 x 16 supermatrices, instead 
of the usual 8x8 for the GOE-GUE crossover @. Here 
we refrain from such a truly heroic enterprise. Instead 
we focus on the regimes of strong and weak spin-orbit 
coupling, and study the intermediate regime by means of 
numerical simulations of the model (^) . 

Before we address the case of strong spin-orbit scatter- 
ing A ^> 1 in the crossover Hamiltonian, we first consider 
the GSE, corresponding to A 2 = AN. In the GSE, the 
wavefunction tf> is a vector of independently Gaussian 
distributed complex numbers. Then, one easily verifies 
that, for large N, the elements of the matrix G of Eq. 
(||) are real random variables, independently distributed, 
with a Gaussian distribution of zero mean and variance 
2/N. Hence G is a random real matrix with distribution g 2 



FIG. 2. Distribution of the orientationally averaged 
gr-factor g 2 — (g 2 + g% + gI)/3 (upper left) and of the ra- 
tios ri2 = I c?i / P2 1 (circles) and r23 = |<?2/<?3| (squares, main 
figure). The solid curves are computed from the theory (fi"l|), 
the data points are numerical simulations of the random ma- 
trix model (||) with N = 200 and A = 7.7. The slight dis- 
crepancy between theory and simulations for ri2 is a finite- iV 
effect; good agreement is obtained with the GSE distribution 
with N — 200 (dotted curve). The lower inset shows {g 2 } 
vs. 1/N for A = 4.3 (diamonds), 6.2 (squares), and 8.1 (open 
circles), together with the theoretical prediction (g 2 ) — 3A~ 2 
for N — » 00 (closed circles). 



The principal 3-factors are the eigenvalues g^ of the prod- 
uct Q = G T G. The distribution of the eigenvalues of such 
a matrix product is known in literature pljfl . It is given 
by Eq. (|) with (g 2 ) = 6/N. 

Let us now turn to the Hamiltonian H(X) for large 
A ^> 1, but still A N 1 / 2 . In that case, spin-rotation 
invariance is broken globally (so that a wavefunction as 
a whole does not have a well-defined spin), but not lo- 
cally; on short length scales, the particle keeps a well- 
defined spin. We then argue that, in the random matrix 
language, one may think of the quaternion wavevector 
ifj as consisting of ~ A 2 ^> 1 components, each with a 
well-defined spin (or "quaternion phase"), but with un- 
correlated spins for each component. The distribution of 
Q is then given by the distribution for the GSE with N 
replaced by a number ~ A 2 |l7j]. We have found that 
the precise correspondence is N — > 2A 2 , by estimating 
the exponential term in the exact distribution, along the 
lines of Ref. [ ^0p^ ] . In order to verify this statement we 
have numerically generated random matrices of the form 
(^). The comparison with the GSE distribution with N 
replaced by 2 A 2 is excellent, see Fig. ||. 

In order to further analyze P(G) for strong spin-orbit 
scattering, we introduce the orientationally averaged 3- 
factor, 



\{gl + gl + gl) = {{25e,/nB\B\) 2 ) i 



(10) 



P{G) oc exp(-A^trG T G'/4). 



(9) 



where the brackets 



indicate an average over all 
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directions of the magnetic field. Further, we introduce 
the ratios ria = I51/52I and r-iz = I52/53I to characterize 
the anisotropy of Q. Changing variables in Eq. (||), we 
find that P(g, ri2, ^23) reads 
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Note that the distribution of ri2 and r 2 3 does not depend 
on (g 2 ) (provided the spin-orbit scattering is sufficiently 
strong). The "g- factor" g z for a magnetic field in the z- 
direction (which is a random direction with respect to the 
principal axes) is given by g z = (Gzz) 1 ^ 2 ■ Its distribution 
follows from Eq. (g) as P(g z ) oc g 2 exp(— 3g 2 /2(g 2 )), in 
agreement with Ref. |J. 

The case of weak spin-orbit scattering can be addressed 
by treating the terms proportional to A in Eq. (|^) as a 
small perturbation. To second order in A we find, 



g = 4 - 4A 2 V 



1 



) 2 ' 
1 J 



(12) 



where A is the mean level spacing and is an antisym- 
metric 3x3 matrix proportional to the matrix elements 
of the perturbation in the eigenbasis of H(0) = S, 

(o/*v)ij = N' 1 / 2 {ip^Auty^euj , where e kij is the anti- 
symmetric tensor. We first consider the change in the 
principal g-factors due to the matrix element cou- 
pling the level e M to a close neighboring level e v where 
v = /i + 1 or /i — 1. (Level repulsion rules out the pos- 
sibility that both levels e M ±i are very close.) In view of 
the energy denominators in Eq. (12), we may expect that 
this contribution is dominant. Taking only the relevant 
matrix element a^ v into account, we find 



.93 = 2, g x = g 2 = 2-\\ 2 {e ll - e v \ 



(13) 



where v = [i ± 1. Since the spacing distribution P(|e M — 
er^l) ~ irA^ 2 ^^ — e„| for small — e„ |t]], we find that 
the distribution P(g) of both g\ and 32 has tails P(g) = 
(3A 2 /2tt)(2 - g) - 2 for 2 - g > A 2 . The main effect of 
contributions from the other energy levels in Eq. ( |l2] ) is 
a reduction of 33 below 2, and a separation of g\ and 
#2. This is illustrated in Fig. |l]. The three regimes of 
weak, intermediate, and strong spin-orbit scattering are 
compared in Fig. |3| using a numerical evaluation of the 
distributions of the three principal g-values. 
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FIG. 3. Distributions of the principal <;-factors gi, g2, gs 
for A = 0.6, 2.0, and 7.7. The data points are obtained from 
numerical simulation of Eq. (M with TV = 100. 



D. C. Ralph, C. T. Black, and M. Tinkham, Phys. Rev. 
Lett. 74, 3241 (1995); ibid 78, 4087 (1997); 
D. Davi dovic and M. Tinkh am, Phys. Rev. Lett. 83, 1644 
(1999); |sond-mat/9910396| . 



D. G. Salinas, S. Gueron, D. C. Ralph, C. T. Black, and 

M. Tinkham, Phys. Rev. B 60, 6137 (1999). 

W. P. Halperin, Rev. Mod. Phys. 58, 533 (1986). 

C. P. Slichter, Principles of Magnetic Resonance 

(Springer, Berlin, 1980). 

K. A. Matveev, L. I. Glazman, and A. I. Larkin, cond- 



mat/0001431 



[7] M. L. Mehta, Random Matrices (Academic, New York, 
1991). 

[8] A. Pandey and M. L. Mehta, Commun. Math. Phys. 87, 
449 (1983). 

[9] J. B. French, V. K. B. Kota, A. Pandey, and S. Tomsovic, 

Ann. Phys. (N. Y.) 181, 198 (1988). 
[10] H.-J. Sommers and S. Iida, Phys. Rev. E 49, 2513 (1994). 
[11] V. I. Fal'ko and K. B. Efetov, Phys. Rev. B 50, 11267 

(1994); Phys. Rev. Lett. 77, 912 (1996). 
[12] S. A. van Langen, P. W. Brouwer, and C. W. J. 

Beenakker, Phys. Rev. E 55, 1 (1997). 
[13] A quaternion is a 2 x 2 matrix q of the form q — 50 1 + 

i y~] . qjffj, where the qj are real numbers (j — 0, 1, 2, 3). 
[14] For example, if the principal axes of Q are labeled ei, £2, 

and £3 = ei x §2, and we apply a static field Be~3, then 

a resonant AC field b oc Re[(g2ei + irje2gi)e~ luJt ], with 

uj — (73 \B/h > 0, will produce spin flips for n = 1 but 

not for r) = —1. 
[15] K. M. Frahm, private communication. 
[16] E. Brezin, S. Hikami, and A. Zee, Nucl. Phys. B 464, 411 

(1996). 

[17] The same is true for the GOE-GUE crossover, where the 
relevant quantity is the "phase-rigidity" p — \ (Tip\ip) | 2 . 
The distribution P(p) for large magnetic fields equals 
P(p) in the GUE, with N replaced by 2a 2 , where a is a 
crossover parameter analogous to A, see Ref. 112]. 



4 



